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We study the statistics of the spacing between Coulomb blockade conductance peaks in quantum dots with 
large dimensionless conductance g. Our starting point is the "universal Hamiltonian" — valid in the g oo 
limit — which includes the charging energy, the single-electron energies (described by random matrix theory), 
and the average exchange interaction. We then calculate the magnitude of the most relevant finite g corrections, 
namely, the effect of surface charge, the "gate" effect, and the fluctuation of the residual e-e interaction. The 
resulting zero-temperature peak spacing distribution has corrections of order A/y^. For typical values of the 
e-e interaction (r^ ~ 1) and simple geometries, theory does indeed predict an asymmetric distribution with a 
significant even/odd effect. The width of the distribution is of order 0.3A, and its dominant feature is a large 
peak for the odd case, reminiscent of the (5-function in the g^ oo limit. We consider finite temperature effects 
next. Only after their inclusion is good agreement with the experimental results obtained. Even relatively low 
temperature causes large modifications in the peak spacing distribution: (a) its peak is dominated by the even 
distribution at AbT ~ 0.3A (at lower T a double peak appears); (b) it becomes more symmetric; (c) the even/odd 
effect is considerably weaker; (d) the 5-function is completely washed-out; and (e) fluctuation of the coupling 
to the leads becomes relevant. Experiments aimed at observing the T = peak spacing distribution should 
therefore be done at T < 0. 1 A /ks for typical values of the e-e interaction. 



PACS numbers: 73.23.Hk, 73.40.Gk, 73.63.Kv 
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I. INTRODUCTION 

The Coulomb Blockade (CB) of electron tunnelin 
of the most studied effects in quantum dots (QDs).E. 
It allows one to probe quantum interference effects in both 
the wavefunction and the energy of interacting electrons. The 
main way in which the latter has been probed is through the 
spacing between adjacent CB conductance peaks. A satisfac- 
tory explanation for the observed CB peak spacing distribu- 
tion (PSD) has, however, remained elusive. In this paper we 
first focus on the T = PSD and find its shape for quantum 
dots containing a few hundred electrons. We then turn to the 
effect of temperature, showing that it is surprisingly large. In 
the end reasonable quantitative agreement between theory and 
experiments is obtained. 

The CB effect occurs when the thermal energy k^T is 
smaller than the charging energy Ec = e^/2C required to 
add an electron to the QD — C is the total capacitance of the 
QD. In that case, electron transport through the QD is blocked 
by energetics, fixing the number of electrons N in the QD. By 
sweeping the voltage Vg of a capacitively coupled gate, this 
CB can be overcome at a particular value Vj^ where the tran- 
sition — > iV + 1 occurs. The conductance G{Vg) shows 
then a series of sharp peaks as a function of Vg as additional 
electrons are added to the QD. At sufficiently low tempera- 
ture, only the ground state (GS) contributes significantly to 
the conductance peak. The position of the CB peak is then 
proportional to the dnange in the GS energy of the QD upon 
adding one electron.cl 

The simplest model used for the description of this phe- 
nomenon assumes a constant e-e interaction — hence the name 
constant interaction (CI) model. The single-particle part of 
the Hamiltonian is described by random matrix theory (RMT) 
under the assumption that the single-particle classical dynam- 



ics is chaotic (or diffusive). As a result, the fluctuation of 
both the conductance peak height and the CB peak spacing 
are given by single-particle RMT statistics. Despite the s 
cess of this model in explaining. the formeilZl — once therma 
and periodic orbit effectsEloEl are iaclu^d .. it f ails drasti- 
cally in describing the observed PSD.HtitiotZy 

Within the CI model, electrons fill the states of the QD in an 
"up-down" scheme due to the spin degeneracy. This implies a 
strong even/odd effect on the PSD. However, none of the ex- 
periments to date have showii-such an effect — though a weak 
even/odd effect was observedll§t3 — suggesting that spin plays 
a more active role. Furthermore, the observed PSD presents a 
Gaussian-like shape (with broader non-Gaussian tails), which 
contradicts the expected Wigner-Dyson distribution from ran- 
dom matrix theory. Finally, the magnitude-pithe width of the 
PSD was questioned. Early experimentsOM found it scaled 
with Ec, which is much bigger than the predicted value, 
^ A, the sinalefmrticle mean level spacing. More recent 
experimentsOtZlLj however, showed that it is indeed of order 
A. 
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experiments, it was suggestedii 
tions were dominated by the e-e interaction itself. Therefore, 
a completely different approach — involving non-perturbative 
methods such as self-consistent Hartree-Fock or exact 
diagonalization — was required. On the other hand, based on 
the fact that a typical QP,'^'|]"tf''I'^|'^ ■V'Vg'j nr'T'^p'' of electrons, 
AT > 1, it was arguedEEjiBmBBlMa that they should 
be described as "good" metals. This implies that the residual 
e-e interactions (i.e., those beyond Ec) are weak and can be 
added to the CI model perturbatively. We shall take the latter 
approach and show that it provides a good description of the 
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experimental data. 

The small parameter in this perturbative approach is 1/g 
with g = E±/ A the dimensionless conductance and Eth the 
Thouless energy (approximately h times the inverse time of 
flight). The condition for the QD to be a good conductor is 
goi\/N '3> l.Ert3 Interaction corrections are classified by their 
order in A /g and successively added to the CI Hamiltonian. 
It then becomes clear why the CI model is wrong: there is a 
zero-order correction — i.e. a correirtkm-of order A — namely, 
the average exchange interaction.EjE^EjtJ Although this is a 
small correction to the total energy of the QD — and so the per- 
turbative approach is justified — it is crucial for properties, like 
the CB peak spacing, that are sensitive to single levels in the 
QD. The zero-order Hamiltonian — hereafter called the. 
stant exchange and interaction (CEI) model — is given byl 



H, 



CEr 



(1) 



where {Eq} are the single-electron energies, Af — CgVg/e de- 
scribes the capacitive coupling to the control gate, Cg is the 
dot-gate capacitance, S is the total spin operator, and J5 is 
the exchange constant. The difference between the CEI and 
CI models is the additional term proportional to 5'^. Because 
it sometimes leads to a GS wi|li,,S|^l, the simple "up-down" 
filling scheme breaks down.EOEZlEl The corresponding PSD 
is completely different from the CI model result (see Ref. ^ 
for a plot). In fact, a GS with S = 1 has been experimen- 
tally observed very recentlyEJ The PSD resulting from the 
CEI model is still, however, in poor agreement with the data. 

Most of the work so far has concentrated on the calculation 
of higher order corrections to the Hamiltonian. The most im- 
portant ones are: (1) the "scrambling" of the spectrum when 
adding an electron to the QD;Ej (2) the fluctuation of the diag- 
onal matrix element of the e-e interaction ;E3 and finally, not 
related to the e-e interaction, (3) the change in the single- 
electron energies when the gate voltage is swept.Ej Surpris- 
ingly, although these corrections have been discussed previ- 
ously in the literature, an explicit calculation of the PSD in- 
cluding all of them has not been done — note however that Ref. 
^ included the first two. Here we present results that include 
all three effects and show that the scrambling and gate effects, 
though domiaaoL, are much smaller than usually assumed in 
the literature.EjEl We also^how that the fluctuation of the off- 
diagonal matrix elementsEjLj introduces a small correction — 
in the regime relevant for the experiments — and can be disre- 
garded in the calculation of the PSD. Despite the substantial 
improvement these corrections introduce, the disagreement 
with the experinieiital results persists. 

Very recently,E!] we pointed out that there is a simple ef- 
fect that has not been taken into account: finite temperature — 
we should mention though that temperatute-effects have been 
discussed in terms of spinless particleslij£3 In our previous 
paperBj we show that in the CEI model the temperature ef- 
fects are more important than in the CI model and that they 
become significant even at k^T ^ O.lA. Since most experi- 
ments were done in the regime k^T ^ 0.3-0.5A — an excep- 
tion is Ref.|8] — our results are crucial for interpreting the ex- 
perimental data. Here, we present those results in more detail 



and extend them. We show that, because of temperature, the 
fluctuation of the coupling to the leads strongly affects the 
PSD. It is not until temperature is introduced that good agree- 
ment with the experimental results is obtained. Furthermore, 
temperature introduces the biggest correction in the regime 
where most experiments were done so far and constitutes the 
main cause of smoothing of the PSD. Consequently, lower 
temperature experiments are required in order to observe the 
actual ground state PSD. 

The rest of the paper is organized as follows: in Section || 
we review the arguments that lead to the CEI model. The lead- 
ing order corrections to this model are introduced in Section 
We calculate the contribution of the off-diagonal terms 
of the e-e interaction and the magnitude of the scrambling ef- 
fect in Section IV and ^ respectively. Numerical results for 
the T = PSD, including all the corrections, are presented in 
Section Jvl We introduce the effect of finit e tem perature in 
Section |VI1|. Finally, we conclude in Section VIII , 



II. THE CEI MODEL 

At low temperature and low bias, only a few energy lev- 
els around the Fermi energy (Ep) are involved in the trans- 
port process. Consequently, an effective Hamiltonian Hqu 
capable of describing the QD in that energy window is all 
that is needed. When considering single-particle properties 
of chaotic (or diffusive) QDs, it is well-known that Hqu can 
be described by random matrix theory (RMT) provided that 
g ^ I. This approach is valid within an energy window 
up to the Thouless energy _Eth- The single-particle Hamilto- 
nian is then "universal", i.e. it only depends on the symmetry 
of the problem — ^broken time-reversal symmetry is assumed 
throughout the paper — and, of course, the energy scale A. 
This approach has proved to be quite successful for under- 
standing the role of mesoscopic fluctuations in transport prop- 
erties of QD (see Refs. ^ and 37 for reviews). 

On the other hand, the treatment of the e-e interaction is 
more subtle. A proper description requires taking into account 
the screening of the bare Coulomb interaction provided by the 
electrons beyond Et^. If the interactions are not so strong, i.e. 
if the gas parameter Vg is small, the screening can be calcu- 
lated using the random phajse approximation. In that case, the 
screened potential readsQcll 



e2 

V:sc(ri,r2) = — + VrF(ri,r2) + F(ri)A + l/(r2)A 



(2) 



where C is the capacitance of the QD, VTF(r) is the Thomas- 
Fermi screened potential and V{r)A is a finite-size screened 
potentialEj — its specific form and origin will be discussed in 
Section ^ Here and throughout this paper we consider two- 
dimensional (2D) quantum dots. The last three terms in Eq. 
(0) are of order A and so much smaller than the first term. 
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This leads to the following Coulomb interaction Hamiltonian 

Him = Ec{n^ -n) + {h-l) ^ C^^aC^^a Xa,f3 



where 



H. 



dridr2*Kri)^';(r2) 

xVrF(n-r-2)*/3(r-2)^'a(n) , (4) 

Xa,f3 = A /drl/(r)5'* (r)^'^(r) and ^^{r) is the eigen- 
function of the single-electron Hamiltonian with eigenvalue 
£a- Fluctuations of the wavefunctions cause both ^ ^ g 
and Xa.p to fluctuate. However, since both matrix elements 
are defined as integrals over the QD's area. A, one could ex- 
pect the main contribution to come from their mean value due 
to self-averaging. It turns out that this is indeed the case and 
that the parameter that controls the relative importance of the 
fluctuations is 1/g (or equivalently l/fcpV^). This suggests 
we can expajidp^ii,, in powers of 1 / g and keep only the dom- 
inant terms .EtJE3 

The zeroth-order term {g oo) in this expansion corre- 
sponds to taking the mean value of the matrix elements, 

where J' = A^'^ J dri dr2 Vtf(t'i - r2)~ A/2 and 



(5) 



(6) 



is the exchange constant, Jo (a;) is a Bessel function (2D case) 
and kp is the Fermi wavevector. We shall see in Section 
that {Xa.p) = A/2(Sq.^. Introducing these mean values 
in Eq. (^ wp-iObtain the "universal" part of the interaction 
HamiltonianEEj 







(0). 



-^Ecfi^-Js 



(7) 



where we have dropped terms linear in n and redefine Ec 
to include all the terms proportional to n^. Here, S — 
J2a,aj^ c^,^'?<t,<t'C„,o.- is the total spin operator of the QD. 
Eq. (0) is the most general Iprm of the interaction Hamil- 
tonian compatible with RMT.EZi Adding the coupling to the 
control gate gives the CEI model Hamiltoruan [Eq. (|l])]. 

Using the explicit expression for Vtf(''')£3 we can write Js 
in terms of the gas parameter Vg, 



5 A Arcsec {rs/^/2) 



(8) 



Note that J5 < A/2 is a fixed quantity. Fluctuations in the 
spectrum of Hcei arise only from {Eq.}. This is a key point 
for understanding its GS: while in the CI model the levels 
are filled in an "up-down" scheme — which leads to a bimodal 
PSD — in the CEI model it is energetically favorable to pro- 
mote an electron to a higher level and gain exchange energy 
whenever the spacing between the top two single-particle lev- 
els is smaller than 2 J5 (N even). This is why the PSD is very 
different from the CI model result. 



III. LEADING ORDER CORRECTIONS TO 
UNIVERSALITY 



In the previous section we considered the simplest model 
for Hqu, which only takes the universal part of the residual 
interactions into account. Here, we include the next order cor- 
rections: (1) the "scrambling" of the spectrum when adding an 
electron to the QD f4 (2) the fluctuation of the diagonal part of 
the e-e interaction ;E3 and (3) the change in.the single-electron 
energies when the gate voltage is swepto Although (2) is a 
correction of order A /g to the Hamiltonian, all three effects 
lead to corrections of order A/y^ to the spacing. We discuss 
now each of them in detail. 



A. Scrambling 

The scrambling effectiil is caused by the presence of the 
potential V{r)Am Eq. (^, which leads to the second term in 
Eq. (^. Its physical origin is quite simple. When an electron 
is added to the QD, the other electrons arrange themselves to 
screen the extra charge. That means that a charge — c/k is pull 
out to the boundaries of the QD. This extra charge creates an 
additional potential, V{r)A, for the electrons inside the QD. 
While in 3-D materials the charge is confined to a small region 
near the surface, in 2-D it is inhomogeneously distributed over 
the whole area of the QD. 

It is straightforward to show that (Xa^p) — VA Sa^p where 
V — Ar^ J drV{r). The mean value thus does not introduce 
any correction to the CEI model — it should be added to Eq- 
The correction comes only from the fluctuations of Xa,p. In 
particular, the main contribulipn (to the spacing) arises from 
Xa,a, with a the top level.ElIu Its variance is given by. 



var(A'„^„) = A2 / dridr2V{n)V{r2){\K{ri)-i'Jr2)h- 



(9) 

withVir)^V{r)-V. 

With corrections of order 1/g included,_the_W|ayefunction 
correlation appearing in Eq. (^ is given byHEJCjCJ 

A'{\mri)^c.{r2)f) = l + k{r,,r2)-kir,)-kir2) 



-k + nB{ri,r2) 



(10) 



with fc(ri, r2) = Jg (fcp |ri - r2|), k{r) = A^'^ J drik{r,ri), 
and k = A^^ j dridr2k{ri,r2)- nB(ri,r2) is a classical 
propagator that contains the contributions of the trajectories 
that reach the boundary of the QD — it is therefore geome- 
try dependent. The latter satisfies / driYlsiri, r2) = with 
i = 1 or 2. The second term in Eq. (|l^ corresponds to 
Berry's result: on scales smaller than the system size ^/A, 
the correlation of chaotic wavefunctions is giving by a ran- 
dom superposition of plane waves.c3 The other terms involv- 
ing k{ri,r2) pipperly account for the normalization of the 
wavefunctions. cJ 
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After substituting (|l^ in we get 



1 1 g but numerically much smaller — so 



A2 /• . . 

var(A'„^a) = \ dridr2V^(ri)V^(r2)fc(ri,r2) 

'dridr2-(/(ri)F(r2)nB(ri,r2). 



(11) 



For quantitative evaluation we consider the case of a ballis- 
tic circplar disc with diffusive boundary conditions. For this 
systemcj 



IIb (t"1,T2) 



1 



- V-1 / rir2 



5=1 



cos(7( 



(12) 

where = r,; (cos 0^, sin 0,^). Thus the last term in (|ll]) yields 
exactly since the potential V(r) is isotropic (in the isolated 
dot case, see Section For a general geometry, the form 
of IIb is not known, and there is no reason a priori to expect 
such a cancellation. We assume that the contribution of this 
term is of the same order as the first and that our final result 
for the var(A'Q,^Q) is correct up to a factor 2.H 

Since fc(ri,r2) ~ l/7rfcp|ri — r2|, a simple dimensional 
analysis of the first term iijii&l- ( pi] ) shows that it is propor- 
tional to lkY\fA so thaOEll 



A2 

vai-(A'Q,Q) = 6oo — • 
9 



(13) 



The same result is valid for var(A'Q.^). Again using the cir- 
cular disc with diffusive boundaries for quantitative estima- 
tion we have 5tii—p fi'^iv-p/R so that g = ^ik-p^/A/2s/TT 
with 7i — 0.38.cJl!] The approximate value of thegeometry- 
dependent coefficient 6oo is calculated in Section \ 



B. Diagonal matrix elements 

Recently, Ref . |9| has shown that the fluctuation of the di- 
agonal terms of iJint leads to a correction to the spacing of the 
same order as the scrambling, despite the fact that the variance 
of the diagonal matrix elements is of order A^/g^ xhe reason 
is that the correction due to these terms involves a sum over 
« g levels. To see this, let us first calculate the variance of the 
diagonal matrix elements. In the zero-range approximation, 
where the short-range screened potential is approximated by 
a (5-function, Vi^'ii')^ &{r)/^A/2, it is given by 



var(AfQ,/3) 



3A2 ln(4/cFVv4) 



47r 



In 9 A 9 
^A^ 

r 



^A^ (15) 



For the double-diagonal matrix element, var(AfQ, Q) — 
4var(Ma,0)£3 Let us now consider the contribution of the 
fluctuation of the diagonal terms to the spacing correspond- 
ing to the transition ^ ^ 0^ i, 



E {Mi 



■+l,/3' 



(16) 



Neglecting, the correlation between the different matrix 
elements£3 we get var(sdiag) — 25var(MQ.;3) cx A^ln g/g.E^ 
The contribution of these terms to the spacing fluctuations is, 
then, of the same order as the scrambling, namely A / y^. 
A similar result can be obtained for all the different spin 
transitions except for ^ ^ 0. In that particular case 
Sdiag — Mn_ jv and therefore var(sdiag) oc /g'^. This is an 
exact result not related to the zero-range approximation for 
Vtf- For this reason, the effect of the fluctuation of the di- 
agonal terms on the (5-function of the PSD is weak, though 
noticeable. 

Diagonal terms produce fluctuations of the energy differ- 
ence between the singlet and triplet states, Eg^Q — Eg^^. 
It is straightforward to verify that var(_Eg^g — -E^^i) — 
2gvaT{Ma^f3). Since the exchange contribution to this energy 
difference is 2 J5, we can think of this as an effective fluctu- 
ation in J5. These fluctuations might lead to a change of the 
ground state in the cases where the triplet and singlet states 
are almost degenerate. 

It is worth commenting that there is a correction of order 
1/(7 to the mean value of the diagonal matrix elements. In the 
zero-range limit it is given by 



'-Cl7^{Sa,sSi3.^ + Sa.',Si3,s) (17) 



with ci = gA^^ J dr [IIb {r,r) — k] cx In ^.0 This correction 
can be included in the definition of Ec and J5 since it has the 
same structure as Eq. (||). For a ballistic disc, this introduces 
a correction of— 0.008A in the latter 

Before closing this subsection we should point out that 
{Ma.p Xa^a) = C(l/.<?^)- In principle this also leads to a 
A /y^ correction to the spacing. Rough estimates suggest that 
it is numerically small, however, and so we neglect it. 



A2 
4^2 



dridr2 



fc(ri-r2)+nB(ri,r2) 



(14) 

where Ma,p = H^pp^^~ {Ha,p,p,a)^°^ and fc(ri-r2) = 
fc(ri — r2) — fc(ri) — fc(r2) + fc. Using^the full expression for 
Vtf leads to similar numerical resultsEj — in that case it is im- 
portant to keep the correlation between the direct and the ex- 
change terms. The dominant contribution in Eq ( p^ comes 
from [fc(ri — r2)]^ — the other terms are of the same order in 



C. "Gate" effect 

The two contributions we have discussed so far are usu- 
ally regarded as intrinsic, in the sense that their origin is 
the e-e interaction of the particles in the QD itself. On the 
other hand, the effect of the gate voltage is usually associ- 
ated with the distortion of the shape of the QD when the gate 
voltage is swept.O Consequently, it appears that it could be 
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independently ijeduced by using, for instance, QDs defined 
lithographicallyllJ or a uniform gate as opposed to a finger- 
shaped one. This is not true however: the gate effect is always 
as big as the scrambling effect. 

The correction to the confinement potential due to a change 
in the gate voltage isQ 



SU{r)^-2Ec5U-A 



(18) 



with SAf = = CI Wile and Wl is the change 

of the electrostatic potential of the i-th gate. The potentials 
V{r) and (r) are related to the solution of the electrostatic 
problem associated with the set of conductors (QD+gates) — 
see next section and Ref. ^ for more details. The former is the 
same potential that appears in Eq. (Q) and causes the scram- 
bling effect. It is clear then that a change in produces the 
same effect as a change of iV. This is reasonable since these 
two effects are opposite faces of the same electrostatic prob- 
lem: a change of the electrostatic potential of the QD produces 
a non-uniform distribution of charge, which creates l/(r)A; 
viceversa, an extra charge must be distributed in the same way 
so that the potential of the QD is uniform. JEq. (|l8]) leads to 
the following correction to the HamiltonianQ 



jj 

gate 



5NX, 



a, 13' 



(19) 



where ^ is defined in terms of y (r) — it is then clear that 
var(A2 fj) = hii^^ I g, with ha a geometry dependent numeri- 
cal coefficient. 

Notice that the change in the shape of the QD can be as- 
sociated with the last term in Eq. (|l8|). For instance, if the 
potentials of two plunger-gates are swept in a way such that 



= SAfi + 5M2 = then the change of the confinement 
potential originates only from that term, and the effect of a 
change in the shape of the QD can be isolated.tj Note, how- 
ever, that this procedure actually tests only the difference be- 
tween V^*^^^ (r) and V'^'^^ (r), which could be smaller that each 
one of them if the two gates are in similar positions with re- 
spect to the QD. 



IV. OFF-DIAGONAL MATRIX ELEMENTS. 

So far we have considered only the contribution of the di- 
agonal terms of H\ai- We found that although the fluctuation 
of each individual matrix element is of the order of A /g, the 
total contribution to the peak spacing is of order A/ y/g. This 
is a consequence of the addition of the ~ g different matrix el- 
ements that contribute to the spacing. Since in principle there 
are many more off-diagonal terms, one might wonder if their 
contribution can also add-up-and result in a significant one that 
should also be included.E3^ We show now that this is not the 
case. 



The first correction to the GS energy due to off-diagonal 
terms appears in second order perturbation theory (the first 
order contribution is zero by definition) 



(20) 



where 1^*5 ) is the GS of the system (described by Hqei) with 
electrons and spin 5", and are the excited states. 

Following Ref. ^ we rewrite the off-diagonal part of the 

Hamiltonian (||) as follows. 



l3>a,i>S 



(21) 



r 



with 



'^a,/37,(5 



H 



(22) 



The sum in Eq. ( pT] ) runs over all configurations in which the 
indices of and c are not fully paired (that is, the configu- 
rations not included in a Hartree-Fock treatment). From this 
form of the Hamiltonian, it is easy to see that conserves 
the total spin.L3 The first term produces only triplet-transitions 
while the second only singlet-transitions. Because of that, the 
states coupled by have the same spin. Then, the energy 



denominator that appears in Eq. ( |20| ) involves differences be- 
tween single-electron energies. The spin rotational invariance 
of H°^^ also implies that the second order correction to the en- 
ergy is the same for all the states in a given spin-multiplet (i.e 
the states with different Sz)- Therefore, we can use the one 
with the maximum value of Sz, the simplest state, throughout 
our calculations. 

(2) 

The difficulty in calculating Eg lies in recognizing which 
terms have to be added coherently, that is to say, which terms 
lead to the same final state This complication arises be- 

cause of: (i) the indices in Eq. (pT|) might be partially paired. 
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so terms within each of the two main terms do not necessar- 
ily produce orthogonal states, or (ii) non-trivial states, like the 
5 = 1 state, can lead to the same final state under the action of 
any of the two main terms in Eq. ( ^1] ) for particular values of 
the indices. In order to avoid the first problem, we explicitly 
take into account all the possible pairing of the indices and 
rewrite Eq. (eTJ) as 



(23) 



with 



1 ^ " t 

+ 2 (."Q,/3,7,Q+^a,/3,7,a) '^a,cr S.^'^/S 



a,0 <y 



(24) 



Hu = 



El t t 

7 

/' t t t t ^ 

i,j>5 

E 9 "a>/3,7,7 4,T4,i('^/3a'^a,T~C;5jC„ j) 



a, 7 



7,/3>a 



(25) 



and JTc as in Eq. (|T]) but with all the indices being different. 
Here, none of the remaining indices are paired and a = —a. 
Notice that Ha contains contributions from both terms in Eq. 

In general, the eigenstates of Hqei are a superposition of 
Hartree-Fock states. However, because S'^ and the single par- 
ticle Hamiltonian commute, they can be classified by their 
occupation numbers For each configuration with Ng 

singly occupied levels, there are 2^= states which have dif- 
ferent values of 5 and S^- It is worth mentioning that for 
Ng > 3, the values of S and Sz are not enough to specify a 
given eigenstate-for instance, for Ng = 3, there are two or- 
thogonal sets of states with S = ^. The determination of the 
precise form of the-spin eigenstates in terms of the HF states 
is not a trivial taskJl3 Nevertheless, for our purpose, it is suffi- 
cient to notice that in order to add coherently two terms must 
lead to final states with the same occupation numbers. With 
this in mind, it is straightforward to check that Ha, Hb and 
He add incoherently. In analyzing each of them however, we 
must proceed in a case by case basis. 

At present, we are not aware of any simple method for eval- 
uating Eq. ( ^0| ) for an arbitrary state l^'^ ). However, it is suf- 
ficient for our purpose to calculate the correction for the cases 
5 = 0, 5 and 1 since those are the most probable values of the 
spin for typical values of the e-e interaction. As an example. 



let us discuss the 5 = case in detail. In this case, l^'^Lo) has 
only doubly occupied levels (up to Ep). It is easy to see that 
the last two terms of Ha give zero when applied on |\E'sLo)- 
The remaining terms of Ha give. 



/3,7 \ a 



■a,/3,7,a ^~ ^a,/3,7 



V2 



(26) 



Note that the sum over a does not affect the final states — this 
is an example of terms that add coherently — while each pair 
{/?, 7} leads to a different final state. The factor \/2 was intro- 
duced to keep the final state properly normalized. Similarly, it 
is easy to check that each of the terms in Hq and Hq lead to 
orthogonal states. Adding up all the contributions, we finally 
get the following expression. 



E 



(2) 

s=o — ■ 

/3,7 
2 



^^a,/3,7,a"l~'^a,/3,7,a 



/3>a,7>(5 ' ' 



2 ' 



(27) 



where a and /3 (7 and S) refer to occupied (empty) levels. 
Notice that the first (coherent) term is absent in Ref. |6[ Fol- 
lowing a similar procedure, we have obtained expressions for 
S — ^ (see Appendix) and S —1 which are too cumbersome 
to be presented here. In those cases, however, the presence 
of single occupied levels introduces some complications, and 
special care must be taken with the terms in which any of 
the indices correspond to one of these levels. It is also im- 
portant to properly take into account the symmetry properties 
of the matrix elements, Ha.p^^.s = Hp^a.s,-,, as it affects the 
variance of Ua^p^^^s and 5. Note than in the particu- 

lar case of the zero-range interaction limit, u^ji^^^s = and 

aQ,/3,7,5 = 2-ffQ,/3,7,5 
(2) 

We calculate Eg numerically in the zero-range limit. We 

f 2) 

first evaluate the correction to the spin gap, A5 = -BgZg — 
f 2) 

EsU- We find that both (As) and rms(A5) are of order 
A /g which, being a higher order correction, can be neglected. 
Nevertheless, we checked that for N — 500, rms(As) is 
~ 4 times smaller than the correction introduced by the di- 
agonal terms. Therefore, we ignored any effect of the off- 
diagonal terms in the occupation of the states. This allows 
us to simply evaluate the correction to the spacing as Soff = 

EnIi,s"+Em-i,s-'^E^n]s- We find that var(soff) oc AVg^ 
so can also be neglected. An explicit evaluation for N = 500 
gives rms(soff) = 0.016A which is much smaller than the fluc- 
tuation introduced by the diagonal terms. The latter result 
corresponds to a generic transition, while for the special case 
^ i -> we get rms(soff) = 0.007A. 
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V. MAGNITUDE OF SCRAMBLING EFFECT 

As we mentioned in the previous section, the origin of the 
potential V{r) in Eq. (||) is the screening charge —e/n. Al- 
though it has been known for a whileEll that this leads to a 
correction of order A/^/g to the Hamiltonian, vw:{Xa,p) = 
booA^/g, a realistic estimate of the magnitude of &oo is still 
lacking. This is particularly important since the effect of the 
scrambling on the PSD goes in the right direction, i.e. it can 
lead to Gaussian-like distribution if it is strong enough. Hejp, 
we show that this has been overestimated in the literatureEfEl 
and that scrambling is not able by itself to explained the ex- 
perimental results. 

The evaluation of 600 for actual geometries is quite difficult. 
The reason is that it involves finding the solution of the eleca 
trostatic field for a set of conductors in a particular geometry Jj 
Following Ref. ^ we write 



where 



Vir) 



An 
8^ 



(28) 



where dz = dz_ — i9z+ (z is the axis perpendicular to the dot), 
I?(ri, is the Green function of the electrostatic problem 
outside the QD including the gates, 

VlV{r,r2)^-Sir~r2), V{r,r2)\reS = (29) 
with S the conducting surfaces and 



An 



dridr2d z,V{ri,r2) 



(30) 



is the total capacitance of the QD. Eq. (^Sj) has a very clear 
physical interpretation if we notice that 



(31) 



is the electrostatic potential outside the QD with (/)(r)|reQD = 
—e/C and (j){r) — over the gate electrodes. Then, 



with 



(32) 



(33) 



the surface charge density in the QD associated with electro- 
static potential (f>ir). Using Eq. (^0|), it is straightforward 
to show that Q — J dra{r) — —c/k; note that this implies 
V = ^. Introducing Eq. (^2|) in Eq. ( pi] ) and using the fact 
that k{ri,r2) ~ 1/Trkp\ri~r2\ we get 



var(A'„,/3) 



dri dro 




(34) 



0'(r-i) = 



\ri - r2\ 



(35) 



is the potential due to the surface charge in the QD. Let us 
now consider the different cases. 



A. Isolated dot 

In this case, the only charge in the system is o'(r) and the 
electrostatic potential (j)'{r) is constant over the QD surface, 
(/>' (r) — <j){r) — — e / C . This allow us to readily obtained 



var(A'c,/3) = 



Ak 



inkpy/A 



C 



(36) 



with a — A^'s J dridr2 |ri — r2|^^. Notice that no particular 
geometry has been assumed so far The first term in the paren- 
thesis can be calculated numerically for arbitrary geometries. 
The second, however, requires the calpilation of the capaci- 
tance. In the case of an ellipsoidal QD,cl 



Ak / b JH 1 
= \\ F(arcsma;, ^) , 




C 



where a (b) is the length of the long (short) axis and m) 
is the elliptic integral of the first kind. For a/b c:^ 1 — 3 we get 
600 — 0.002. It is clear then that Tms{Xa,f3) — 0.04A/^ is 
smaller than usually assumed. In fact, for N = 500, this value 
corresponds to Sx — y/boo/g = 0.018 in the parametric ap- 
proach to the scrambling, about a factor of 10 less than taken 
inRefs. ||and|3|. 



B. Dot with gates 

Since experiments are certainly done in the presence of 
gate electrodes, a careful calculation should take them into 
account. Looking at Eq. (|6|), it is tempting to simply re- 
place C by its experimental value. However, the above cal- 
culation is only valid for an isolated QD, as we explicitly as- 
sumed that (T(r) was the only charge in the systems. For a 
real QD, the induced charge on the gates has to be consid- 
ered. Then, (f){r) — (f)' (r) + 4>" (r) where 4>"{r) is the po- 
tential created by the induced charge on the gates. Defining 
<f" = A-^ J dr0"(r) and (3 = (Qcf")-^ J drcr(r)0"(r) we 
find 



A2 



An 



Airkpy/A 



C 



l + (2- 



(38) 

Notice that the value of the ratio / C can now be obtained 
from the experimental data since C is the capacitance for the 
actual geometry. 

Since (3^1 and cj)" > (the sign of the induced charge 
is the opposite of Q), it is evident that using the isolated dot 
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result [Eq. ([36[)] gives an upper limit to var(A'c j3) when eval- 
uated with the experimental parameters. An estimate of is 
obtained as follows. Let Qi be the charge of the z-th gate, then 
(f)" ~ Qi/diK, with di the distance between the centers of 
charge of the QD and the i-th gate. Since Qi =—Cg{—e/C), 
where C* is the dot-«-th gate capacitance, it turns out that 



r / (e/ C)^J:^ CI /d,K. Then, 



^0.5 



(39) 



C e/C 

where we used that typicall}0 X^i ^'g/C^ 0.5 and denote by 
d the average distance between the center of the dot and the 
gates. For the data of Ref. |l^ we estimate 5oo ^0.005 with an 
upper limit of 0.01; thus, even with gate effects-included, our 
estimate is smaller than values used previously .Eic^ 



VI. GROUND STATE PEAK SPACING DISTRIBUTION 

We now use numerics for the evaluation of the PSD. At T = 
(temperature effects will be discussed in the next section) the 
peak spacing is given by 

= [eS+' (M' ) - E^s i^')]- [E^s W E^s' m (40) 

where EQg{JV) is the GS energy of the QD with N electrons 
excluding the ch^ing energy term and J\f is the correspond- 
ing gate voltageEd Including all the leading order corrections 
to the CEI model, the Hamiltonian readsO 



^QD — H, 



CEr 



diag 

4 E 



+ E <.^P,Ain~M)X^,p + SMXl^] (41) 

where the variance of p ^ g and ^ are given by Eq. ( [l5| ) 
and Eq. (|8]), respectively. Their mean values are included in 
the definition of Ec and J5, so that {H^ p j s) — i^i p) — 
0. The second term in Eq. ( |4l| ) includes only the diagonal 
terms of the residual interaction, and 5Af is taken with respect 
to some fixed state, for instance the state with N electrons. 
The GS energies were obtained by minimizing the energy with 
respect to the occupation numbers for A^ — 1, N, 1, and 
N+2 (and the corresponding A/^, Af', andAf"). We only kept 
two consecutive spacings for each realization of the single- 
particle Hamiltonian. It is worth mentioning again that Eq. 
( |4l| ) is defined in a window up to Ei^, so that only g levels 

were considered. 

The parameters we use are: LLlg — 0.38^/ N/2 (which cor- 
responds to a disc geometry) ;EE1I (2) the upper limit for the 
value of boo — 0.01, and (3) 611 = boo- Only transitions with 
1(55*1 = I were taken into account since in the absence of spin- 
orbit interaction (or if it is small), the transitions where the 
change in the GS spin is bigger than i, appear as "missing" 
peaks in the conductance and are not included in the experi- 
ment. 




0.75 



0.5 I 

spacing [A] 

3 




0.5 I 1.5 2 -0.5 
spacing [A] 



0.5 
spacing [A] 



FIG. 1: Top panel: Ground state peak spacing distribution obtained 
from Eq. (0) with = 500 {g ^ 6) and Js = 0.28A (r^ = 1). The 
solid (dashed) line correspond to TV odd (even). The origin of the 
horizontal axis corresponds to sm ~ Js /2. The CEI model result is 
included in the inset for comparison. Note that all the sharp features, 
including the 5-function, are washed-out but the long tail for large 
spacing persists in the even case. Bottom panel: Comparison of the 
even (left)and odd (right) distributions for N = 200 (solid line) and 
— 1000 (dashed line). Clearly, residual effects are stronger for N 
odd. 



Numerical results for A^ = 200, 500, and 1000 are shown 
in Fig. |l| The horizontal axis corresponds to sn — Js/2, so 
that the origin agrees with the Hartree-Fock result of Ref. ^ 
Corrections to the CEI model clearly smear-out all the sharp 
features of the PSD. The former (5-function (see inset in Fig. 
|l]) is now a finite peak but still constitutes the dominant feature 



TABLE I: Comparison of the 1 / ^ corrections for different number 
of electrons, A^.rWe use the relation g ~ 0.27v'iV, which is valid for 
a disc geometry Ej 



iV 


g 




rms(A/a,/3) 


rms(sjv) 


200 


3.8 


0.051A 


0.031A 


0.313A 


500 


6.1 


0.041A 


0.020A 


0.308A 


1000 


8.6 


0.034A 


0.015A 


0.305A 



9 



of the PSD. As expected, the additional corrections increase 
the rm.s. of the spacing with respect to its value in the CEI 
model (0.28A) (see Table |). Notice that 1/^ corrections 
mainly affect the odd distribution — the even distribution re- 
mains essentially unaltered. This is evident in the comparison 
between the N — 200 and N — 1000 cases (see bottom panel in 
Fig |l]). Therefore, both the scrambling and the gate effect are 
the dominant effects as they are the ones that most affect the 
(5-function — the effect of the diagonal terms on t he (5-f unction 
is higher order in 1 /g (see discussion in Section [II B ). 

It is important to emphasize that besides the smearing 
caused by the additional corrections, there is still a noticeable 
even/odd effect. This theory then does indeed predict such an 
effect at low T. To what extent this is valid at higher T is 
discussed next. 



VII. TEMPERATURE EFFECTS 

So far we have ignored thermal excitations and calculated 
the PSD in terms of GS energies (Fig.|l]). This will remain a 
good approximation so long as k^T-^S, where 6 is the energy 
difference between the GS and the first excited state. In that 
case, the contribution from the excited states can be ignored. 
This has been an implicit assumption in most previous work 
(note however Refs. 



16 and 35 ). We will show now that in the 



CEI model the condition for achieving fceT^C (5 is much more 
restrictive than in the CI model and that in fact temperature 
effects are crucial for understanding the experimental data.Ej 
There are two simple reasons to expect stronger temper- 
ature effects in the presence of exchange. First, there is a 
change in the occupation of the excited states. For example, 
assuming N even, we have S — \ Ae— 2 J5 1, with Ae the single- 
particle energy spacing between the two top levels. A signifi- 
cant occupation of the first excited level occurs when k^T^d. 
Since we are assuming k^T ^ A, this implies (5 ^ A. While 
in the CI model the probability for that to occur is small due to 
level repulsion, this is not so in the presence oiexchange (i.e. 
Ae ~ 2Js is much more likely than Ae ~ 0).Ej Second, at fi- 
nite temperature the peak position involves the change in free 
energy of the QD upon adding ajparticle. Then, as we show 
below, the entropy contributiorocZcj leads to a shift of the 
peak position that depends on the spin transition in the QD. 



given by 



Wo 



Here, 



(43) 

(1) is the equilibrium probability that the QD con- 



tains N electrons, (2) iJQol^'f ) = Cjl^f ) so that "j" labels 
the many-body states of the QD, (3) Feq(j|iV) is the condi- 
tional probability that the eigenstate j is occupied given that 
the QD contains N electrons, and (4) /(e) = {1+ cxp[(e — 
EF)/kBT]}^^. Since near the peak onlx the states with A^— 1 
and N electrons are relevant, we havecZl 



N 



exp (— 



N 



exp (— f2Ar)+cxp (— r^Tv-i) 



(44) 



with N the grand-canonical potential and the canonical 
free energy of the QD. To make the dependence on Af explicit, 
let us denote by {Ej} the eigenenergies of Hqo without the 
charging energy term and define SJ^ = (A^ — 4) —J^. Then, 



E. 



N-l 



2Ec5N 



N-\ 



■-N-l 



E. 



N 



Ef-^ + kBTln 



^eq(jliV) 



F,^{i\N~l) 



{Q.EcSU 



(45) 



with Ztv the canonical partition function. Note that Eq. ( p3| ) 
is valid for any i and j. The contribution of the transition 
i io the conductance reaches its maximum when f{J-N — 
JFjv-i)[l — /(ej — 6;)] peaks, namely when 



Ef = E. 



N 



'E, 



k^T 



In 



^eq(jliV) 

F,S\N~l) 



-2Ec5M. (46) 



B. Ground state dominated transitions 

In the particular case where the transition between GS dom- 
inates, and taking the spin degeneracy into account, the CB 
peak position is given by 



A. General approach 

We now proceed with a detailed calculation. Let us con- 
sider the regime F ^ k^T, A ^ Ec, where F is the total 
width of a level in the QD. Near the CB peak corresponding 
to — 1 ^ A^ transition, the linear conductance is given 



^ pN 



F^F 



Wa 



(42) 



with Fq'^' the partial width of the single-particle level a due 
to tunneling to the left (right) lead and Wa a weight factor 



Eri = Et 



N 

GS' 



k^T 



In 



1 



25, 



N-l 
GS 



1 



-2Ec5N. (47) 



We see that the peak is shifted with respect to its position at 
T — Q bAL^M|-amount depending on the change of the spin of 
the QD.calica Except for a factor ^ in front of the entropic 
term, Eq. (^) corresponds to replacing Eq^ by Tn in the 
usual condition for the peak position, which is what we would 
naively expect at finite temperature. Because the rm.s. of the 
PSD is ^ 0.3A (see Fig. 1), this shift is significant even for 
k^T ^ O.lA and cannot be neglected. Notice that we have 
not made any assumptions about the Hamiltonian of the QD 
so far — except that close to the conductance peak it depends 
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on the gate voltage only through the charging term.0 While 
in the CI model this introduces-only a constant shift between 
the even and odd distributions^ in the CEl model it changes 
the shape of both distributions since different spin transitions 
contribute to each one. Also, one should note that this en- 
tropic effect shifts the energy E^s in the same direction as the 
exchange interaction. Then, we should expect an effect on the 
PSD similar to the one corresponding to an effective increase 
of Js- 



C. Peak conductance 

It is important to point out that the magnitude of the 
on-peak conductance is renormalized because of the spin 
degeneracy OcZl The reason is that Feq{j\N), the overlap 



Ar-1\ 



and the value of x [1— /(e^— e^)] at its 
maximum depend on the particular spin transition involved. 
In the simplest case when only the GS is relevant, we get 



F,^{j\N){P^^[l-f{e,~c,)]}U = 



and 



1 



E 

S',k':S^,k;> 



(48) 



25' + 1 ifn„ = 



25+1 ifria^l 

(49) 

where S' (S), S'^ {Sz) and k' (fc) are the quaotum numbers 
associated to the state with N (N—l) particles. H At low tem- 
perature, most transitions correspond to the first case in Eq. 
(^) when S'>S and to the second when S' < S. Using that, 
we finally get 



G, 



peak " 



r^r 



2e2 



nksT VMS' 



with 



2(5' + 5) + 3 



4 (V25' + 1 + ^/25+1)^ 



(50) 



(51) 



Then, the average conductance peak depends not only on the 
average coupling to the leads but also on the probability of the 
transition S' — i.e., it depends on Js and on the statistics 
of the single-particle spectrum and so on magnetic field. 

This is relevant for a quantitative understanding of the low- 
temperature ^havior of a = 1- (Gpeak)GOE/(Gpeak)GUE in 

closed QDs.lII At T = 0, since higher spin is more likely in 
the GOE case and since A is smaller the bigger the spins in- 
volved, this renormalization leads to values of a larger than 
0.25 — how much larger, of course, depends on Js- At finite 
temperature, when several transitions contribute to the con- 
ductance, it might also lead to values larger than 0.25 and 
could explain the small deviation observed at low temperature 



(ksT<QAA) in Ref. Notice that either the CI i 
dephasing processes lead to values smaller than 0.25.1 



Ae,-2J 

1 s 



Ae,+Ae.-2J 

1 2 s 



Ae, 



Ae,+Ae,-3J 

1 2 s 



Ae, 



FIG. 2: Schematic representation of some low-energy states consid- 
ered in the calculation of the PSD. The scheme at the top (bottom) 
corresponds to A'^ even (odd). Below each state is shown its energy 
difference from the left most one within the CEI model. Aei and 
Ae2 are the first and second level spacing respectively. 



D. Several state case 

As we mentioned above, in the general case more than one 
transition contributes to the conductance, and the CB peak 
position must be determined by maximizing Eq.(42) with re- 
spect to Af. For arbitrary T, this requires the calculation of 
all possible transitions between the eigenstates of Hqu with 
N~l and N electrons. For simplicity, we restrict ourselves to 
low enouah temperature so that only a few excited states are 
relevanto Therefore, we kept 6 states in the even case and 4 
in the odd one. We checked numerically that at k^T — 0.3A 
and for Js = 0.28A, the occupation of these states is, on av- 
erage, 99.4% (98.3%) for N even (odd), being smaller than 
98% (92%) only about 1% of the time. In any case, the effect 
of temperature can only be underestimated since, in general, 
different transition leads to different spacing which in turns 
leads to an smearing of the PSD. Figure ^ shows some of the 
energy states considered in the calculation of the PSD. For 
even, the lowest states with 5 = and 5 = 1 are the dominant 
states with an occupation of 52.4% and 39.6% respectively 
at fefiT = 0.3A and Js = 0.28A. In the odd case, those are 
the lowest states with 5 = i and 5 = |, with 80% and 7.5% 



respectively. 

The upper panel in Fig. H shows the CEI model PSD for 
non-zero temperature. Besides the expected smearing of the 
sharp features and their shift due to the entropic term in (|47|), 
there are two important new effects: (a) Temperature alone is 
able to completely wash-out the J-function, making the odd 
distribution broader. Note in addition the long tail for large 
spacings; we show below that the latter is not simply thermal 
broadening, (b) The even distribution develops a peak at small 
spacings — in particular, the maximum of the total distribution 
is dominated by the even distribution, in sharp contrast to what 
occurs at r = 0. This strongly reduces the relative weight of 
the long tail in the even case, and the distribution becomes less 
asymmetric. Actually, the long tail is only slightly affected by 
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FIG. 3: Top panel: Finite temperature CB peak spacing distribution 
corresponding to tlie CEI model [Eq. (|l|)] witii Js = 0.28A (r^ = 1) 
and kaT = 0.3 A. Notice that the (5-function in the odd distribution 
(solid line) is completely smeared out by temperature and that the 
even distribution (dashed line) develops a peak. Bottom panel: Same 
but neglecting the fluctuation of the coupling to the leads. Note that 
the long tail in the odd distribution (solid line) for large spacing is 
absent, and that the width of the peak in the even distribution is sig- 
nificantly reduced. 



ent levels (Fq,). This should be particularly important when 
the GS and the first excited state are almost degenerate. The 
bottom panel of Fig. ^ shows the PSD assuming F^ = cte. and 
using the same parameters as before. It is evident that much of 
the broadening observed in the top panel is not directly caused 
by temperature but by the fluctuation of Fq. 

One of the most important differences is the absence of the 
long tail for large spacing in the odd distribution. This can be 
easily understood as follows. First, let us note that the sharp 
jump in the T = PSD at Js results from the transitions in- 
volving 5=i: ^ i^O, 0^ 1, 1^ 5^0 
and 1 ^ -i ^ 1. The spacing in each case is Js, Ae2 — Js, 
Aei—Js and Aei+Ae2— 3 Js, respectively. It is easy to show 
that at T = 0, the conditions on Aei and Ae2 for each transi- 
tion to occur, that is for the GS to have the appropriate spin, 
leads to a spacing < Js in all the cases. This result is a conse- 
quence of the "yes-no" conditions required at T = 0. At T 7^ 0, 
those conditions are relaxed and the last three transitions can 
lead to a spacing bigger than Js- This thermal broadening 
is responsible of the disappearance of the (5-function. How- 
ever, because the reahzations contributing to that part of the 
distribution have Asi ~ 2 J5, the thermal factors of the dif- 
ferent transitions are very similar to each other. Consequently, 
the relative strength of the couplings can overcome them: the 
peak position is dominated by the most strongly coupled level, 
which might coiTespond to the larger spacing. This explains 
the larger tail observed when the fluctuation of F^ is taken 
into account. 

Similarly, the width of the peak of the even distribution is 
strongly affected. This clearly indicates that fluctuations of 
the wavefuntions of the QD strongly modify the PSD. 

So far we have discussed temperature effects in the context 
of the CEI model. It is surprising that, even at this level of 
approximation, only a weak even/odd effect or asymmetry is 
expected for A:bT>0.3A. 



temperature as it corresponds to large values of the single- 
particle spacing. 

The peak in the even distribution arises from cases where 
S=l and S = states are (almost) degenerate. It corresponds 
to the sharp discontinuity at the origin in the T = PSD — 
where both spin states significantly contribute to the conduc- 
tance. Then, all the transitions with Aei ~ 2 J5 collapse into a 
single (average) value for the spacing, which leads to a peak in 
the PSD. According to Eq. (^6|) the corresponding CB peaks 
are shifted by ^ zt^k^T ln(4/2), which gives a total shift of 
kjiT In 2 for the peak in the PSD. 

Note that the rm.s. of the distribution is reduced by temper- 
ature. In fact, we found that it decreases monotonically from 
T = 0. 

The fact that more than one transition contributes to the 
conductance implies that the peak position also depends on 
the relative strength of the coupling to the leads of the differ- 



E. Corrections to CEI model 

Results including the leading order corrections are shown 
in Fig. Q for fcB^^ O.lA and 0.3A with the same parameters 
as in Fig. |l]. The additional fluctuations increase the broaden- 
ing of the distribution. At low temperature, the effects of the 
residual interactions are dominant — though the T-induced- 
peakin the even distribution is evident. For fcBr = 0.3A, how- 
ever, temperature is the main effect (compare with Fig. |]). In 
this case, the even/odd effect is weaker but stifl noticeable — it 
should be kept in mind that the experimental noise may con- 
ttibute significantly to weaken this effect. Alfo, notice that 
the PSD is not Gaussian. A detailed analysisll3 of the experi- 
mental data of Ref. [l^ shows that this is indeed the case. In 
fact, the agreement between these data and the PSD shown 
in the lower panel of Fig. ^ is good, both quaUtatively and 
quantitatively. 
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FIG. 4: Finite temperature CB peak spacing distribution correspond- 
ing to the Hamiltonian Eq. (0) with Js = 0.28A, TV = 500, and 
boo = 0.01. The top (bottom) panel corresponds to ksT — O.IA 
(0.3A) and the solid (dashed) line to A'' odd (even) 



VIII. CONCLUSIONS 

In this work we have calculated the Coulomb blockade peak 
spacing distribution including the most representative leading 
order corrections (up to A/ y/g) to the CEl model as well as 
the effect of finite temperature. 

At T = 0, our results show that the PSD still presents a 
clear signature of an even/odd effect. Even though it is much 
weaker than the effect predicted by the CI model, it is defi- 
nitely big enough to be observable. No sharp features remain, 
and the peak in the odd distribution (the former (5-function) is 
still the dominant characteristic. Also, the distribution is far 
from being Gaussian and its width is ^ 0.3A. This number, 
however, depends on the value of Js and on many geometry- 
dependent parameters that could vary a bit for the actual QD. 
On the other hand, it could be argued that the RPA approach 
used in the calculation of the screening of the Coulomb poten- 
tial is not appropriate for ~ 1. We think, however, that the 
essential ingredients are captured by this approach and that 
any correction that would arise from a more accurate calcu- 
lation could be included by a renormalization of Ec and Js- 
Such a renormalization could have an important impact both 



in the shape of the distribution and in the strength of eveao/odd 
effect since they are quite sensitive to the value of J^.a An 
experimental determination of Js should, then, be a high pri- 
ority. 

For r 7^ the picture is quite different. At k^T 0.3A, 
the roles of the even and odd distribution are inverted: it is 
the even distribution that shows a peak, while in the odd one 
the (5-function is washed-out. This is very important, since 
the absence of the (5-function has been one of the puzzles in 
interpreting the experimental results. The final distribution 
is closer to a Gaussian-like shape, and the even/odd effect is 
much weaker. Here, the main effect in the distribution comes 
from the temperature. One important consequence of the finite 
temperature is that the fluctuation of the coupling to the leads 
becomes relevant and substantially contributes to the broad- 
ening of the PSD. Both the shape and the rm.s. of the dis- 
tribution agree with the data in Ref. |l5| We should mention 
however, that we fail in reproducing the long tail for small 
spacing. This could be due to (a) relevance of higher excited 
states — note that this part of the distribution corresponds to 
cases where the single-electron spacing is very small — and 
the consequent effect of the fluctuation of the couplings, or 
(b) mixing of the top levels caused by the off-diagonal terms 
of the interaction — here our second order perturbation the- 
ory fails — or by the off-diagonal terms of the scrambling. At 
lower T both distributions show a peak, which is a clear ob- 
servable feature. 

It is important to point out that our results do not explain the 
data of either Ref. |7| or |l8|. Nevertheless, in both cases there 
are some elements to think that this is not a "failure" of the 
model. In Ref. Hthe interpretation of the transport process 
itself is not clearH For example, the width of the CB peak is 
not controlled by temperature though its shape corresponds to 
a thermally assisted process. In Ref. |l8| temperature effects 
are negligible due to both k^iT ^ 0.05A and Js ^ 0.25A. 
However, the single-particle dynamics is not fully chaotic be- 
cause of the regular shape of the QD. Therefore, the effect of 
regular orbits and of the presence of regular islands in phase 
space, must be considered. This could enhance the contribu- 
tion of IIb in Eqs. ( pi] ) and (jl^) and lead to larger fluctuations 
of both Ma,p and Xa,p- This subject is quite complex and we 
leave it for future work. Nevertheless it is important to men- 
tion that the rm.s. in Ref. |l8| is of order 0.4A and that there 
is a weak even/ odd effect. 
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APPENDIX A: THE DIFFUSIVE CASE 

Similar results for the leading order corrections can be 
obtained in the case of diffusive QDs after the appropriate 
change of some definitions.Q First, the dimensionless conduc- 
tance is now given by g = hji / A where 71 is the smallest 
non-zero eigenvalue of the diffusion equation, — I?V^/„(r) = 
7n fn{i'), supplemented with von Neumann boundary condi- 
tions. D ~ is the diffusion constant and i the mean 
free path. Second, fc(ri,r2) = exp(— r/i?) Jq (fcpr) with r ~ 
|ri — r2|. Third, the propagator nB(ri, r2) is replaced by its 
diffusive counterpart 



(Al) 



In order to obtain a numerical value for the fluctuation of the 
different matrix elements, a specific geometry must be as- 
sumed. For a disc of radius R we get. 



R? ' 



J-m iPm.n) — 



and 



with 



fra,n (^) -^ri 



COS mq 
sin md 



Ar,-. 



f3„ 



(l + <5™,o)^ i?y^/32„ „-m2|J,„(/3„,„)| ' 



(A2) 



(A3) 



(A4) 



Here, Jm{x) is a Bessel function and J^i^) its derivative. 
Note that the first non-zero eigenvalue, 71 1 is proportional to 
the square of the first zero of J{, /3i 1 = 1.84. This means that 
the last mode to relax is the one with the smallest number of 
nodes in the angular direction and none in the radial direction 
(except for the origin), /i,i(r) cx coscj) Jf{Pi^ir / R). Then, 
g = nDf3lJR'^A cx C/RVn. For N - 500 and £ - R/2 this 
gives gss 13. 



1. Scrambling 

The main contribution in this case comes from the second 
term in Eq. dill). Because there is no reason to assume a 1/r 



decay of the wavefunction correlation in a general case, we 
cannot use the same approach wejused in Section Instead, 
we assume a disc geometry to getfl 



^00 = 7 



71,1 



}_ V- sin^(/3o,„) 
0.004 



drl^(r)/o^„(r) 



(A5) 



which is very close to the value obtained for the ballistic case. 
Here we used V(r) = l/4[l-(r/i?)2]i Thus MW^X^fi) re- 
mains of the same order as in the ballistic case. 

2. Diagonal elements 

In this case, the terms in Eq. (^4|) that involve k{r\ — Jfa) 
are small (assuming 1<^\J~A) and can be neglected. TherBEJ 

var(A/,,;5) = ^ y"dridr2 [nD(ri,r2)]'. (A6) 

Assuming a disc geometry we get ^wijxla.ji) — c^I^ j'^g^ 
with 



C2 



7i,i 



00 00 



EE 



0.13 



(A7) 



Once again, the numerical result is similar to the one we ob- 
tained for the ballistic case. For instance, using N ^ 500 and 
t^Rjl, we getrms(A/a,;3) W0.014A, which should be com- 
pared with the second row in Table |. 



APPENDIX B: SECOND ORDER CORRECTION FOR S = i 

We give here an explicit expression for the second order 
correction to the energy due to the off-diagonal terms for the 
case of 5*= i. 
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E' 



(2) 



\ ^ 3 |Ma,/3,7,l|^+|aa,/3,7,l|^ , 5a,P . ^ 



Sa5^/3('^'"a,/3,7,a + aa,/3,7,a)(l 



3 I |2 
+ 4 |'"1,/3,1,7 + CH,/3,1,7| 



/3>Q,7 



/3,7 



2(e/3-£7 



|Ec«^/3(3Ma,/3,l,a + aa,/3,l,a) + a/3,i,i,i + ai,i,i,/3| |Ea(3Wa,l,7,a + aa,l,7,a)| 



16(£;3-£l) 



4(£i-£^ 



(Bl) 



where a and /3 (7 and 5) correspond to doubly occupied 
(empty) levels of the S = ^ state — except for the term with 
E' in which case the singly occupied level labeled "1" is in- 



cluded. Note that there are several terms which add coher- 
ently. 
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One should note, however, that the covariance between the matrix 
elements Ma^ji and is of order A^/(j'^ (see Ref. ^ for de- 

tails) which, in general, also leads to a correction of order A / ^ 
to the spacing. For a diffusive QD, this term is small due to the ex- 
ponential decay of fc(ri — r2), but it might be important in cases 
with strong semiclassical effects. 

With this definition, sjv has the classical result 2Ec subtracted. 
That is, Sjv = (Ae) for the odd (even) case in the CI model. 
This is no longer true once the gate effect is introduced. How- 
ever, taking this into account would introduce a correction of or- 
der (A/Ec) X Tms{Xa^i3), which is negligible. 
The quantum numbers S and are not enough to specify a many- 
body state if the number of singly occupied levels is bigger than 
2, even for a given set {wq } of occupation numbers (see Ref. p3| ). 
Notice, however, that a always first decreases from its T = 
value and then increases again. This is related to the fact that the 
contribution from the excited states come firts in the GOE case. 
Because the fluctuations of the single particle levels, this is of 
course always violated for some (very) rare realizations. This pre- 
vents us from describing the asymptotic regime on the left side of 
the PSD in Fig. D and p 



